Simulation of partial entanglement with non-signaling resources 
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With the goal of gaining a deeper understanding of quantum non-locality, we decompose quantum correlations 
into more elementary non-local correlations. We show that the correlations of all pure entangled states of two 
qubits can be simulated without communication, hence using only non-signaling resources. Our simulation 
model works in two steps. First, we decompose the quantum correlations into a local and a non-local part. 
Second, we present a model for simulating the nonlocal part using only non-signaling resources. In our model 
partially entangled states require more nonlocal resources than maximally entangled states, but the less the state 
is entangled, the less frequently must the nonlocal resources be used. 

PACS numbers: 



I. INTRODUCTION 

Quantum correlations are very peculiar, especially those vi- 
olating some Bell inequality jit] . Gaining a deeper insight into 
such nonlocal quantum correlation is a grand challenge. Chil- 
dren gain understanding of how their toys function by disman- 
tling them into pieces. In the present paper we follow a similar 
approach by decomposing the quantum correlations into sim- 
pler, more elementary, nonlocal correlations. 

This work is part of the general research program that looks 
for nonlocal models compatible and incompatible with quan- 
tum predictions. The goal is to find out what is essential in 
quantum correlations. Note that we do not claim that Nature 
functions as our model. Nevertheless we believe that finding 
the minimal resources sufficient to simulate quantum corre- 
lations, and studying the computational power that they offer 
10, S B Htl, provide enUghtening insights into the quantum 
world. 

In the last years, two different ways of decomposing quan- 
tum correlations have been proposed. The first one, due to 
Elitzur, Popescu and Rohrlich (whence EPR-2 !) |@], consists 
in decomposing some quantum correlations into a local and a 
non-local part. A second approach consists in the simulation 
of entanglement with the help of some non-local resource, e.g. 
classical communication or a non-Iocal box. While communi- 
cation models |0,llt] give insight to quantum correlations from 
the point of view of communication complexity, we believe 
that models using only no-signaling resources are more 
relevant from a physical point of view, since it is most unlikely 
that Nature uses any form of communication ifioll . In this pa- 
per, we shall combine for the first time both approaches, and 
prove that all pure entangled states of two qubits can be simu- 
lated using only no-signaling resources, i.e. without commu- 
nication. 

The approach we follow works in two steps. First, in the 
EPR-2 spirit, we decompose the quantum correlation Pq cor- 
responding to von Neumann measurements performed on pure 
entangled states of two qubits = cos 6*100) + sin 6*111) 

into a statistical mixture of a local correlation Pl and a non- 
local correlation P^l |6(|: 



Pq=Pl{0)Pl + {l-pLmPNL ■ (1) 

The weight PL (6*) is thus a measure of the locality of the state 
|V'(6')). In particular, for any maximally entangled state of 
two-qubits one has pl{9 ^ tt/A) = la], a result that holds 
true for maximally entangled states in any dimension ifTHl . 
Note that in general the probability distribution P^ l does not 
need to be quantum, but is restricted to no-signaling correla- 
tions by construction. 

Then, we provide a simulation of the nonlocal correlation 
Pnl using only nonlocal, but non-signaling resources. Ac- 
cordingly, in order to simulate Pq, it suffice to simulate Pl 
with probability pl{0), which requires only shared random- 
ness (but no nonlocal resources), and to simulate Pnl with 
the complementary probability 1—^^(6*). As expected, the 
less the quantum state IV'(^)) is entangled, the smaller the 
weight pl{Q) of the local correlation Q [l2ll . Consequently, 
the simulation of a less entangled state requires less frequent 
use of nonlocal resources; in particular for separable states 
Pl{S = 0) = 1. However, not much is known about the 
nonlocal resources needed to simulate the nonlocal part of 
quantum correlations, i.e. to simulate pnl- Reference ||^ 
presented a simulation of the quantum correlation for the spe- 
cial case of maximally entangled qubit pairs (9 = tt/A) using 
only one nonlocal box, the so-called PR-box ifTsll . For non- 
maximally entangled qubit states, very few is known. To our 
knowledge, the only known result is that one PR-box is not 
sufficient for simulating slightly entangled states lITill . This 
result shows that entanglement and non-locality are different 
resources, as also suggested by other works lITsl [l6l [l7i [Tsll . 

In this paper we use a decomposition of the form ([l]), re- 
cently presented in Ref. lfT2ll . which is optimal under some 
general assumption, and present a simulation of the corre- 
sponding nonlocal correlation Pnl for arbitrarily entangled 
two qubit states. This simulation requires finitely many non- 
local boxes, though no claim of optimality can be made. For 
pedagogical reasons, the paper is organized as follows. After 
introducing the general framework in Section 2, and briefly 
reviewing the case of maximal entanglement in Section 3, we 



2 



present in Section 4 a preliminary model for simulating par- 
tially entangled qubit states, without using any decomposition 
into local and nonlocal parts. This allows us to introduce the 
two main ingredients of our model: first the technique of cor- 
related local flips; second the Millionaire box, a generalization 
of tiie PR-box. Then in Section 5, we briefly recall the de- 
composition into local and non-local parts presented in lfT2ll . 
and explain how our preliminary simulation model can be ex- 
tended to simulate the nonlocal part P/v l of the model of Ref. 
m. Finafly we give some conclusions and perspectives. 

II. GENERAL FRAMEWORK 

Formally, a correlation is a conditional probability distri- 
bution P{aj3\ah), where a, (3 denote the outcomes observed 
by Alice and Bob when they perform measurements labeled 
by a and h. Here, measurements are conveniently repre- 
sented as vectors on the Bloch sphere, since we focus on 
von Neumann measurements on qubits. A correlation is non- 
signalling if and only if Alice and Bob's marginals Ma and 
Mb, are independent of the partner's input: Ma does not de- 
pend on b and Mb does not depend on a. For binary outcomes 
(a, (3 G { — 1, +1}), the correlations are conveniently written 
as 

P{a, (i\a, 6) i (l + aMA{a) + /3Mb(6) + a(iC{d, 5)) (2) 
where 

Ma (a) = ^aP{a,(3\d,h) 

a,l3 

MB{b) = Y,pP{a,(3\a,h) , (3) 

are the local marginals, and 

C(a, b)^Yl af3Pia, f3\a, b) (4) 

a, 13 

is the correlation term. Here we shall focus on pure entan- 
gled states of two qubits \ip{d)) = cos0|OO) + sin 6*111}, 
9 e]0,7r/4]. Thus the quantum correlation PQ{a(3\ab) is 
given by 

MA{d)=ca, , MB{b)^cb, 
C{d, b) = a^bz + s{a.j,bx - Uyby) , (5) 

where c = cos 29 and s = sin 29. 

Now, we would like to decompose the correlations Pq into 
simpler ones, such that 

PQ{a,l3\a,b)^ f dXPxia,(]\a,b) , (6) 



where d\ is a normalized measure. In a simulation model, 
two ingredients are required: first, a non-local resources for 
creating the elementary correlations P\; second, a strategy 
(represented by the A's) for judiciously combining them. In 
this paper, we provide such a decomposition. The remarkable 
feature of our model is that the elementary nonlocal correla- 
tions are obtained without communication, that is using only 
no-signaling resources. 

III. MAXIMALLY ENTANGLED STATE 

Let us briefly the simple case of maximal entanglement, 
i.e. 9 = 7r/4. In this case the marginals vanish, MA{d) = 
Mb (b) = 0, and the correlation takes the simple scalar prod- 
uct form C{d,b) ~ a ■ b ll32ll . In reference a model sim- 
ulating this correlations is presented. This model uses as re- 
sources only shared randomness and one PR-box, which sat- 
isfies the relation a (B b = xy, where x, y are Alice's and 
Bob's input bits, and a, b their outcome bits. In general non- 
local boxes provide some elementary nonlocal correlations. 
They are elementary in that they allow only for a limited (usu- 
ally finite) number of inputs and outputs and they are ex- 
tremal points in the convex set of nonsignaling correlations 
ifigll . They are nonlocal in the sense that they violate some 
Bell inequality. Importantly, they do not allow signalling, that 
is the statistics of the local outcomes (i.e. the marginals) are 
independent from the other parties inputs. This model demon- 
strates that the resource needed to simulate maximally entan- 
gled qubit pairs is surprisingly simple. Indeed, what could be 
simpler than a © 5 = xyl 

IV. PARTIALLY ENTANGLED STATES: PRELIMINARY 
MODEL 

We now turn to partially entangled states of two qubits. 
In general the marginals MA{d) and MB{b) do not vanish. 
However, in the case of two parties and binary outcomes, it 
is proven that all extremal nonlocal boxes have vanishing (or 
deterministic) marginals ll20l I21I1 . This explains in part why 
it is difficult to simulate partially entangled states. In order 
to circumvent this difficulty we introduce now the concept of 
correlated local flips. 

A. Correlated local flips. 

Let us consider an arbitrary probability distribution 

Po{a,p\d,b) = \{^ + cc[3Co{d,b)) . (7) 

with vanishing marginals and correlation term Co (a, b). Now, 
Alice and Bob perform local flips on the probability distribu- 
tion Pq, that is, Alice (Bob) flips her (his) output —1 with a 
probability fa (fb), while the output +1 is left untouched. Af- 
ter this processing, also called a Z-channel, the marginals are 
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clearly biased towards +1. Let us now assume that /t > fa 
and that the flips of Alice and Bob are both determined by a 
shared random variable A uniformly distributed in [0, 1]. Al- 
ice and Bob flip their -1 outcome if and only if A < /a and 
^ < fb, respectively. The resulting probability distribution 
reads 

Pfia,l3\a,b)= (8) 

i(l + afa + 13 fb + al3{fa + (1 - fb)Co{S, b))) . 

It should be pointed out that the flips fa and must be corre- 
lated; this will be crucial in the following. Note also that ev- 
ery probabiUty distribution P{a, (3) = \ + uMa + PM b + 
a(3C) with Mb > Ma can be generated in this way. 

B. Preliminary model, step 1 

We just described a technique for creating a probability dis- 
tribution Pf with nontrivial (i.e. non vanishing) marginals, 
starting from an initial probability distribution Pq which had 
trivial marginals. Now the intuition is the following: since 
correlation with trivial marginals seem to be easier to create 
with standard nonlocal resource (such as PR-boxes), let us do 
the identification Pf = Pq and find out what is the required 
initial distribution Pq. For partially entangled states of two- 
qubits (Pq given by Q), this leads to 



C. The Millionaire box. 

Two millionaires challenge each other: who is richer ? 
Since millionaires are in general quite reluctant to reveal how 
much money they own, they prefer to use the Millionaire-box 
(M-box) 1I23II . a nonlocal two-input two-output non-local box. 
The two outputs a,b are binary, (a, b G {0, 1}), and are locally 
random in order to ensure no-signaling. The two inputs x, y 
can be chosen in the continuous interval [0, 1]. The M-box is 
characterized by the following relation: 

a®b=[x <y] , (12) 

where [X] denotes the logical value of X: [X] = when 
X is true. Note that the M-box admits an infinite number 
of possible inputs. So, both millionaires input the amount of 
money they own x,y into the machine; the parity of the out- 
puts (a © b) indicates the winner. Fortunately, the M-box is 
also useful to physicists, as will be shown in the next section. 
Note that the M-box is a generalization of the PR-box; in case 
the inputs x, y are binary, the M-box is simply equivalent to a 
PR-box (given here by x{y + 1) = a ® & © 1). It is also worth 
mentioning that the M-box reaches the no-signaling bound of 
all the Bell inequalities Inn22 introduced in fl^. An inter- 
esting question is whether all (bipartite) non-local boxes with 
two-outcomes ll20ll2lll can be simulated with one M-box. In- 
deed, a detailed study of the non-local properties of the M-box 
would be relevant, but is beyond the scope of this paper. 



fa = caz , fb = cbz , Co = a - B (9) 

where 

B = isb,,-sby,b,~c)/il-cb,) . (10) 

Note that ||i3|| = 1. Remarkably, B corresponds to Bob's 
original measurement setting b moved one step back on the 
Hardy ladder iB. 

Consequently the problem of simulating correlations orig- 
inating from von Neumann measurements on partially entan- 
gled states reduces to the problem of simulating the unbiased 
probability distribution 

Pa = ^{l + al3a ■ B) . (11) 

Such a "scalar product" correlation can be reproduced with a 
single bit of communication ^ or with a single PR-box j^. 
However, there is a caveat: Alice and Bob must know wether 
bz > o-z (as assumed above) or if on the contrary > 6^ ! 
This is due to the fact that the local flips must be correlated. 
Note that in the case > b^, the initial probability distribu- 
tion is given by Pq = j{l + a(3A ■ b), where A is defined 
similarly to equation (fTol i. 

At first sight it may seem that a resource solving this prob- 
lem will lead to signaling, because it would reveal a relation- 
ship between Alice's and Bob's measurements. Remarkably, 
this is not the case. Next, we show that a no-signaling (non- 
local) resource known as the MilUonaire box is exactly the 
tool we need. 



D. Preliminary model, step 2 

As shown above, the technique of local flips allows one to 
recover the correlation of partially entangled states, under the 
condition that bz > dz (or Oz > bz)- But how do Alice and 
Bob know whether bz > az or az > bz 1 The M-box can 
overcome this problem. 

Alice and Bob share two PR-boxes for creating "scalar 
product" correlations (see Fig. [T]i; from now on we call these 
CGMP-boxes 19|]. The first one is used to create the correla- 
tion given by the scalar product a - B, i.e. corresponding to the 
case bz > az and the second one for the scalar product A ■ b, 
i.e. for the case > bz- Local flips are then performed. At 
this point, Alice and Bob have each got two possible outputs 
ai , a2 and /3i , /?2, but don't know which one to use, since they 
don't know whether < b^ or b^ < az- 

Next, they input the z-component of their measurement set- 
ting (respectively and b^ 1330) into the M-box, and get out- 
puts a and 5. It is clear that, for the simulation to succeed, 
the final output of Alice and Bob, a and f3, should be equal to 
Q!i,/3i if az > bz, and equal to a2,/32 if bz > az- Mathemati- 
cally this translates into the following expression 



a © /3 = (a © b){ai © © (a © 6 © l)(a2 © /?2) ■ (13) 
Developing the previous equation, one gets 
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FIG. 1: Preliminary model. Simulating partial entanglement with- 
out communication. The model requires four PR-boxes and a 
Millionaire-box (M-box). The first two PR-boxes create "scalar 
product" correlations (CGMP-boxes). Then the M-box "selects" the 
correct CGMP-box, without revealing any relation between Alice's 
and Bob's measurement settings (i.e. without signaling). Finally, two 
additional PR-boxes are required for computing the correct outputs. 



a©/? = a(ai©a2)©Q;2®6(/3i©/?2)®/32 

©a(/3i © /32) © © 02) , (14) 

which contains some local terms, as well as some non-local 
terms. Remarkably, the non-local terms (second line of equa- 
tion ( fT4l i) are simply obtained by using two supplementary 
PR-boxes, 03 © 63 = a{(3i © (32), and 04 © 64 = h{ai © 
(see Fig. [T]i. 

So finally, using four PR-boxes (two CGMP-boxes and 
two additional PR-boxes) and one M-box, one can simulate 
the correlation of any partially entangled state of two qubits. 
Whether the M-box can be replaced by a finite number of PR- 
boxes (or more generally with a nonlocal box having a finite 
number of possible inputs) is an interesting open question. 
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FIG. 2: Sketch of the main model. 1) Alice and Bob inside the slice 
(see text); both CGMP-boxes can be used. 2) Alice inside. Bob out- 
side; then indeed Uz < fez. 3) Alice outside. Bob inside. 4) Alice 
and Bob outside. Depending wether < or > b^, the M-box 
selects the correct CGMP-box. Note that we have omitted the two 
additional PR-boxes (PR3, PR4 of Fig. 2). 



V. PARTIALLY ENTANGLED STATES : MAIN MODEL, 
INTEGRATING EPR-2 

We are now ready to present our model, combining the pre- 
liminary model (presented in the previous section) and the de- 
composition of Ref |fl2|l, into local and non-local parts (i.e. of 
the form ([TJ). The decomposition is the following: 

Pl{9) = 1-s 

Pl = i(l + a/(a,))(l + /3/(5,)) (15) 

Pnl = ^{l + aF{a,) + PF{b,) + apG{dh)) 

where f{x) = sgn(a;)min(l, j^lxj) , F{x) ~ 7(^2^ ~ 

(1 - s)f{x)) , and G{ab) = Uxb^ - a„h,, + ^[a^bz - (1 - 
s)f{az)f{bz)] ■ We refer the reader to 1I12II for further details. 

Let us point out two important features of decomposition 
(fTsT l First, the weight of the local part pl{S) — 1 — s is a 



monotonic decreasing function of 6, i.e. of the degree of en- 
tanglement of the state \il){0)). Note also that Pi (6') = 1 — sis 
optimal under the assumption that Pl depends only on az and 
bz- Second, the non-local part Pnl depends on the measure- 
ment settings. More precisely, when the measurement setting 
of Alice is such that < (1 — s)/c (i.e. inside a slice of 
the Bloch sphere around the equator), her local marginal van- 
ishes; and similarly for Bob. On the contrary, when the mea- 
surement setting lies outside the slice, the marginal is biased. 
When both the settings of Alice and Bob are found inside the 
slice, the correlation reduces to a simple scalar product with 
trivial marginals. 

The simulation of P/vl is very similar to that presented 
above, thus we only describe Alice's and Bob's strategies. As 
previously, the required non-local resources are two CGMP- 
boxes, an M-box and two additional PR-boxes. After estab- 
lishing non-local correlations with both CGMP-boxes, Alice 
and Bob perform local flips. Finally they use two additional 
PR-boxes to compute the correct output (see Fig. 

Alice proceeds as follows. When her setting is inside the 
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slice (qz < (1 — s)/c), she inputs according to a into both 
CGMP-boxes, and does not perform any local flip (fa = 0). 
When her setting is outside the slice, she inputs the first 
CGMP-box according to A = {sQx, sa.y, c — az)/(l — Uzc) 
and the second CGMP-box according to a. Then she biases 
her output towards outcome +1 with probability fa ~ F{az)- 
Bob proceeds almost similarly. When his measurement set- 
ting is inside the slice (b^ < (1 — s) / c), he inputs both CGMP- 
boxes according to b' = {bx, —by, —bz). When his setting is 
outside the slice, he inputs the first CGMP-box according to 
b' and the second according to i? = {sbx, —sby, 62 — c)/(l — 
bzc). Then he biases his output with probability fb = F(bz)- 



VI. CONCLUSION AND OUTLOOK 

By dismantling the quantum correlations of partially entan- 
gled states of two-qubits into more elementary nonlocal but 
no-signaling correlations, we gained insight into the quantum 
world. We showed that the correlations of all pure entangled 
states of two qubits (under von Neumann measurements) can 
be simulated using only non-signaling resources, hence with- 
out communication. Our decomposition is likely not to be 
optimal in the sense that there might exist more economical 
models. Still, there are already two lessons we learn from 
the present decomposition. First, the less the quantum state 
is entangled, the less frequently one needs to use nonlocal 
resources to simulate it; as intuition suggests. Next, when- 
ever one needs nonlocal resources, then these are definitively 
larger for (at least some) partially entangled states than for the 
maximally entangled state; indeed this is proven for slightly 
entangled states [14:], but is still an open question for close to 
maximally entangled states. Hence, in counting the resources 



required to simulate two-qubit states, one should distinguish 
between the required amount of nonlocal resources and the 
frequency at which one has to use them. 

It is interesting to establish the following connection with 
Leggett's approach to quantum correlation f^?], which re- 
cently attracted quite some attention Ii26i , i27. .28, 29, 30, 31il. 
In models a la Leggett one assumes that the elementary cor- 
relations, contrary to PR-boxes, have nontrivial marginals; 
Leggett's original idea is that each qubit, when analysed in- 
dividually, appears to be always in a pure state, see ll25ll29ll . 
However, one can prove that any such model, with elemen- 
tary correlation having nontrivial marginals, fails to reproduce 
the quantum correlation of maximally entangled states of two- 
qubits ll29i[3lll . This is a kind of converse to the present paper 
in which we show that it is especially hard to simulated at the 
same time nonlocal correlations and non-vanishing marginals. 

Among the open questions, we like to underline the follow- 
ing one. How could one prove that a decomposition is mini- 
mal? As said, this question has two sides. Minimality of the 
resources, and minimality of the frequency at which one has 
to use them. Our experience suggests that the first aspect is 
an especially difficult problem. The second aspect looks more 
promising: it seems natural to conjecture that an EPR2-type 
decomposition with pNL{d) = 1 — cos 26 should exist 1I12I1 . 
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